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The electronic structure and equilibrium geometry of La2/3Sri/3Mn03 are studied theoretically by means 
of density functional calculations. The doping is treated by introducing holes and a compensating jellium 
background. The results for the local density approximation (LDA) agree with previous LDA calculations, 
with an equilibrium volume 5.3 % too small and with both majority and minority spin states present at the 
Fermi level for the relaxed system. The generalised gradient approximation (GGA) offers a qualitatively 
improved description of the system, with a more realistic volume, and a half-metallic behaviour for the 
relaxed structure, which enables studies needing theoretical relaxation. The ideal Mn02 -terminated (001) 
surface is then described with explicit doping. 
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1 Introduction Ferromagnetic perovskite manganites attract much attention because of interesting 
properties like the colossal magnetoresistance found in some phases ' and the high degree of spin po- 
larization of the charge carriers in other phases. The optimally doped case of La2/3Sri/3Mn03 (LSMO) 
reaches complete polarisation at the Fermi level, the so-called half-metallic state. ^ These features make 
manganites very appealing for spin-polarized current injection in new devices for spintronics, such as tun- 
nel junctions. First-principles calculations could be very helpful in the understanding and characterisation 
of such devices. 

Descriptions of the electronic structure of bulk LaMnOa based compounds within density-functional 
theory have already been documented, ^ mainly based on LDA. Electrons can be highly coiTelated in man- 
ganites, and it is known that the electronic characteristics of some of their phases would demand a much 
more explicit consideration of the non-local exchange (Hund's rule in the d-shell of the Mn ions) than 
what provided by present-day GGAs. LSMO, however, is metallic and thus much more amenable to a 
Kohn-Sham description, as already shown by Pickett. ^ 

Two important limitations have been found in the LDA description of LSMO. (i) The theoretical vol- 
ume obtained from LDA is underestimated by 5.3 %, the material losing its half-metallic character in its 
relaxed geometry. " This limitation has been used to argue for the need of a treatment of correlation effects 
beyond LDA, such as the LDAh-U approach,"*'^ even if the point is still debated.^ The flawed volume in 
itself does not prevent realistic studies of the bulk of the pure material, since the geometry can be taken 
directly from experiment. It becomes a very important limitation, however, for any more complex situation 
(surface, interface, defect), for which the first-principles relaxation of the structure is needed before any 
other property can be analysed, (ii) The Sr doping is hard to describe, demanding large unit cells and sta- 
tistical sampling of configurations, or alternative alloy descriptions like the coherent-phase approximation, 
which enormously complicate any study. 

In this work we explore the GGA description of LSMO, observing that a much better estimate of the 
geometry is obtained, which maintains a half-metallic character It is important to stress that this is not 
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a fundamental result, since we do not know whether this character is retained for the right reason. This 
result is important for practical purposes, since it allows the use of GGA DFT for structural optimisations 
of LSMO and its surfaces, interfaces, defects etc. We have also tested the performance of a homogeneous 
approximation to hole doping, replacing Sr substitutions by an homogeneous background of compensating 
charge, as recently tried in other perovskite oxides.^ Finally, the GGA description is tried on the ideal 
(001) surface of LSMO. The homogeneous background would have no sense for this calculation since the 
Sr substitution should be distributed only within the solid and not in vacuum. Hence, we choose simple 
reaUsations of explicit doping for this latter study. 

2 Method We have performed density functional theory (DFT) calculations using the SIESTA method 
which is based on pseudopotentials and numerical localized atomic orbitals as basis sets. The calculations 
were performed using two approximations to DFT, namely, the local density approximation (LDA) and 
the generalized gradient approximation (GGA), both in the spin-polarised version. For LDA we use 
the Ceperley- Alder exchange-correlation potential, whereas for the GGA calculations, we use the the 
Perdew-Burke-Ernzerhof (PBE) scheme. ^■^ 

Core electrons were replaced by norm-conserving pseudopotentials using the scalar-relativistic TroulUer- 
Martins scheme^* in the Kleinman-Bylander factorised form, with nonlinear core corrections. La and 
Mn are dehcate species for pseudopotentials. The 5s^ 5p^ 5cf 4/" and 3s^3p^ 3d^ 4/° configurations 
were used as reference for La and Mn respectively. Notice that the 5s^ and 5p^ semicore states of La and 
the corresponding 3s^ and 3p^ states of Mn are explicitly included in the calculations. The core radii for 
the s, p, d and / components of the pseudopotential are 1.85 ao, 2.20 ao, 3.10 ao, 1.40 ao for La, with a 
radius for the core correction (rpc) of 1.50 ao. The radii for Mn are L40 ao, 1.90 ao, 1.50 ao, and L90 
ao with rpc=0.70 ao. For O and Sr we use the pseudopotentials from previous work. For the GGA O 
pseudopotential we use a slight partial core correction to avoid the GGA instability close to the nucleus. 

The basis set consists of strictly localised atomic orbitals at the double-^ polarised level. The semicore 
states are treated at the single-(^ level. The radial shape of the orbitals is determined variationally under 
a confining pressure of 0.2 GPa, using cubic bulk LaMnOs as reference for La and Mn, and using a H2O 
molecule for O. Pseudopotentials and the basis were tested by comparing our results for cubic LaMnOs 
with those in the hterature both for LDA and GGA^-''. 

The Sr doping is treated in two ways. Firstly, by adding the desired concentration of holes to the pristine 
system, plus a compensating homogeneous negative background.^ This is justified by the main finding that 
the doping produces mainly a change in band filling with little electronic or geometric distortion otherwise. 
Secondly, the results for that approximation are compared with calculations with expUcit Sr doping. 

3 Half-metallicity in bulk LSMO In this section we consider the effect of the density functional in 
the lattice constant and the half-metallic character of bulk LSMO. Although LDA is able to reproduce the 
half-metallic behaviour for the experimental pseudo-cubic lattice parameter^' (ciexp =3.89 A), this is not 
the case when the theoretical equilibrium lattice constant is used (a^^"^ =3.82 A), as can be observed 
in Fig.l (a). The origin of the problem is the underestimation of the lattice parameter %), for which 
the minority-spin band gap is smaller. Within GGA, however, the half-metallic character is maintained 
for both the experimental and the theoretical lattice constants (a'^^^ =3.91 A), which display a much 
better agreement (an overestimation of ~ 0.5%) and thus allow the use of GGA for structure optimisation. 
The dependence of the half-metaUicity on volume is further explored in Fig. 1 (b), where it is shown that 
for values of the lattice constant below ~ 3.80 A the half-metallic character is lost. GGA is used in the 
calculations presented in the rest of the paper. 

The GGA band structure of LSMO is shown in Fig. 2. The bands near the Fermi level are mainly due 
to hybridization of the oxygen 2p and the manganese 3d orbitals. The conduction bands at the Fermi level 
correspond to the occupied Cg majority spin states. The peaks in the DOS at ~ 1.5 eV below Ep and ^ 2 
eV above Ep correspond to the t2g states for the majority (occupied) and minority (unoccupied) spins. 
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Fig. 1 (a) Comparison of the density of electronic 
states for the LDA and GGA approximations. Ma- 
jority (minority) spin states are depicted in the up- 
per (lower) graph of each panel. Solid lines corre- 
spond to the theoretical equilibrium lattice constant, 
while dashed lines are for the pseudo-cubic experi- 
mental value, (b) Density of electron states around 
the Fermi level for three different values of the lat- 
tice constant, within the GGA approximation. 




4 Surface properties For the study of the (001) surface of LSMO we have chosen Mn02 termina- 
tion. A slab with 6 unit cells in the perpendicular direction (12 layers of Mn02-La(Sr)0 planes) and 4x2 
unit cells in the parallel directions was used to model the system. The chosen thickness ensures that the 
bulk properties are recovered at the center of the slab, and the in-plane repetition was introduced for the 
description of explicit 1/3 doping, by substituting one La by Sr in every three. 
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Fig. 3 Projected Density of States for different 
Mn02 layers in the (001) slab. Majority (minor- 
ity) spin is shown in the upper (lower) panel. The 
shifts in the tgg peaks as well as the partial occu- 
pation of the minority spin for the surface layer are 
highlighted by arrows. 
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In Fig. 3 we plot the projected density of states (PDOS) on the Mn02 layers at different depths. The 
PDOS on a layer is obtained by summing the PDOS for all the basis orbitals associated to that layer, and 

PDOS^iE) = C;,,C,nS,f,5{E - En), 

71, W 

where E^ and V'n('^) ~ C'/^n'/'^l'^) are the Kohn-Sham eigenvalues and eigenvectors, respectively, 
the latter expanded on the atomic basis {4>^i{r)}, with overlap S^^i- The surface layer introduces a strong 
peak near the Fermi level mainly for the majority spin but with some weight on the minority spin as well 
that breaks the half-metallicity of the system. This deviation is very localised close to the surface. For 
the subsurface Mn02 layer the PDOS already resembles that of the bulk and for the central layer the half- 
metallic character is fully recovered. This surface contamination of the spin polarisation was found in other 
calculations of doped manganites and is related to magnetic reconstructions at the surface due to changes 
in the electronic occupation for the and dxy states. The t2g states in the surface are shifted upward 
by about 1 eV and overlap with the Cg states at the Fermi level. This changes the magnetic interaction and, 
although the system remains ferromagnetic, the magnetic moments near the surface are slightly reduced 
i^'iAfjiB as compared to '^3.5/iB in bulk). This reduction could be related to the experimental evidence 
of a reduction in the magnetization at the surface of LSMO, but is in contradiction with similar theoretical 
calculations for doped manganites that predict an increase in the magnetic moment at the surface. A full 
characterisation of the surface state and its dependence with termination and surface relaxation will be 
presented elsewhere. 



5 Conclusions We have studied the bulk properties of hole-doped La2/3Sri/3Mn03 comparing LDA 
and GGA (PBE). LDA gives a poor description of both the structure (lattice parameter) and the electronic 
properties (half-metallicity) of the system. GGA reproduces the half-metallic character and gives lattice 
constants in good agreement with experiment. The fact that the full (variable cell) relaxation of the bulk un- 
der GGA retains the half-metallic character allows obtaining relaxed geometries within GGA for surfaces, 
interfaces, or nanostructures of any sort. 
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